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Abstract. We classify recurrent configurations of the sandpile model on the complete bipartite 
graph K m<n in which one designated vertex is a sink. We present a bijection from these recurrent 
configurations to decorated parallelogram polyominoes whose bounding box is a m x n rectangle. 
Several special types of recurrent configurations and their properties via this bijection are examined. 
For example, recurrent configurations whose sum of heights is minimal are shown to correspond to 
polyominoes of least area. Two other classes of recurrent configurations are shown to be related to 
bicomposition matrices, a matrix analogue of set partitions, and (2+2)-free partially ordered sets. 

A canonical toppling process for recurrent configurations gives rise to a path within the associated 
parallelogram polyominoes. This path bounces off the external edges of the polyomino, and is 
reminiscent of Haglund's well-known bounce statistic for Dyck paths. We define a collection of 
polynomials that we call q, t-Narayana polynomials, defined to be the generating function of the 
bistatistic (area, parabounce) on the set of parallelogram polyominoes, akin to the (area, hagbounce) 
bistatistic defined on Dyck paths in Haglund |18j . In doing so, we have extended a bistatistic of 
Egge, Haglund, Kremer and Killpatrick 12 to the set of parallelogram polyominoes. This is one 
answer to their question concerning extensions to other combinatorial objects. 

We conjecture the q, t-Narayana polynomials to be symmetric and prove this conjecture for 
numerous special cases. We also show a relationship between Haglunds (area, hagbounce) statistic 
on Dyck paths, and our bistatistic (area, parabounce) on a sub-collection of those parallelogram 
polyominoes living in a (n + 1) x n rectangle. 



1. Introduction 

The abelian sandpile model [9] is a discrete diffusion model whose states are distributions of 
grains on the vertices of a general directed graph. A vertex of a graph is stable if the number of 
grains at the vertex is strictly smaller than its out-degree, otherwise it is called unstable. In addition 
to a randomised addition of grains, the dynamics of this model requires that an unstable vertex 
donate a grain to each of its direct successors. This is the so-called toppling or avalanche process of 
the model. This process defines a Markov chain on stable states. In the particular case of graphs, 
there exist many bijections from recurrent states to spanning trees of the same graph [91 [61 HI Q] . 
Another rich result concerning the model shows that the recurrent states, together with a certain 
binary operation, form an abelian group called the sandpile group (see e.g. [8]). 

In this paper we classify recurrent configurations of the sandpile model on the complete bipartite 
graph K m>n in which a single designated vertex is a sink. Cori & Poulalhon [7j classified recurrent 
configurations of the complete {k + l)-partite graph K Vx ,__ Pk i where the sink is the single vertex 
that is connected to all other vertices. Their paper showed that such recurrent configurations 
can be related to a generalisation of parking functions, and also exhibited a rich connection to a 
Lukasiewicz language which aided in the enumeration of their new generalised parking functions. 

This paper complements their paper by breaking the symmetry of K Ptqt \ and showing how recur- 
rent configurations of K mn that have one sink can be interpreted as parallelogram (or staircase) 
polyominoes. We study several different types of recurrent configurations, which are illustrated in 
Figure [2j We show that minimal configurations (those whose sum of heights is minimal) correspond 
to ribbon parallelogram polyominoes (those having minimal area). Almost non-zero configurations 
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are configurations in which one distinguished vertex is allowed to be empty, and all others are 
non-empty. 

Bicomposition matrices are defined to be square matrices whose non-empty entries partition the 
set {1, . . . , n} and for which there are no rows or columns consisting of only empty sets. We show 
that minimal almost non-zero ('minanz') configurations on K n ^ n are in one-to-one correspondence 
with bicomposition matrices. A special class of these minanz recurrent configurations are in one-to- 
one correspondence with upper triangular bicomposition matrices. Such matrices were the subject 
of several recent papers that presented surprising connections between five seemingly disparate 
combinatorial objects: (2+2)-free posets, Stoimenow matchings, ascent sequences, permutations 
avoiding a length-3 bivincular pattern, and a class of upper triangular matrices [31 [11] . We present 
the composition of this new correspondence and the correspondence between matrices and posets to 
show how the heights of a configuration can be read from the corresponding poset. It was through 
these structures that we first noticed a connection to the sandpile model and the initial motivation 
behind this paper. 

During the last two decades, a series of papers have examined a power series C n (q,t) that has 
become known as the q, i-Catalan function (or polynomial). This power series was introduced by 
Garsia and Haiman [15] and has important links to algebraic geometry and representation theory. 
In the original paper they showed that two special cases of this polynomial had combinatorial 
significance. The first was that C n (q, 1) was the generating function of the area statistic over Dyck 

paths having semi-length n. The second was that q^ 2 ^C n (q, 1/q) is the nth q-Catalan number. 

Haglund introduced a new statistic 'bounce' (we will call this hagbounce so as not to confuse it 
with another statistic) and conjectured that C n (q,t) was the generating function of the bistatistic 
(area, hagbounce) on the set of all Dyck path of semi- length n. Garsia and Haglund [TH1 Q3] proved 
this conjecture using methods from the theory of symmetric functions. (See Haglund [T7] for a 
concise overview of these results.) Egge, Haglund, Kremer & Killpatrick |12] asked if the lattice path 
statistics for C n (q,t) can be extended, in a way which preserves the rich combinatorial structure, 
to related combinatorial objects. 

We present a pair of statistics (area, parabounce) on parallelogram polyominoes which is one 
answer to their question. We call the resulting polynomials q, i-Narayana polynomials since they 
specialise to the Narayana numbers for the case q = t = 1. The bivariate generating function for 
this pair of statistics F m ^ n (q,t) appears to be symmetric in both q,t and m,n. We conjecture and 
discuss the symmetry of this bistatistic and prove it for numerous special cases. 

The outline of the paper is as follows: Section 2 attends to the classification of recurrent con- 
figurations and a proves a bijection from this set to a set of decorated parallelogram polyominoes. 
Section 3 looks at several different types of (minimal) recurrent configurations and their relationship 
to other structures. The different collections we look at are illustrated and summarised in Figure [2] 
In Section 4 we define the q, t-Narayana polynomials and investigate their symmetry. First we will 
define parallelogram polyominoes and some aspects thereof. 

A parallelogram polyomino is a polyomino such that the intersection with every line of slope 
37r/4 is a connected segment. Parallelogram polyominoes may be described in several ways and 
in this paper we will see three distinct ways of doing this with each serving a different purpose. 
Let us write Para mjn for the collection of parallelogram polyominoes whose bounding rectangle is 
[0, m] x [0, n}. These polyominoes have been studied in several papers, see for example [20] . 
Every polyomino V £ Para m>n can be uniquely described by a pair (upper("P), lower("P)) of paths 
which begin at (0, 0), end at (m, n), take steps in the set {n = (0, 1), e = (1, 0)}, and only touch at 
their endpoints. 

Given V G Para mjn , let area('P) be the area of V. Let us define uarea('P) and larea('P) to be the 
areas of the regions above and below V in its bounding m x n rectangle, so larea("P) + area("P) + 
uarea("P) = ran. Let Ribbon mjn be the set of polyominoes V in Para miT1 whose area is m + n — 1, 
the minimal value, which we will call ribbon polyominoes. 
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Figure 1. Three examples. 

Example 1.1. The first polyomino V\ in Figure^is a parallelogram poly omino. It can be described 
by the pair of paths (upper(T'i), lower('Pi)) = (neneeneenee, eeeneneenen). The second polyomino V2 
is not a parallelogram polyomino, as is witnessed by the two disjoint segments the dashed line of 
slope 37r/4 meets. Polyomino V3 is a ribbon (parallelogram) polyomino. 

Given V G Para m)n , let Bounce('P) be the unique path from (m — l,n) to (0,0) which is defined 
as follows: Starting from (m — l,n), the path goes south until it encounters a vertex on the lower 
path of V . The path then turns to the west and continues straight until it encounters a vertex on 
the upper path of V . The path turns again to the south until it encounters a vertex on the lower 
path of V, and so on, until it reaches (0,0). 

Let bounce(T') = (ci,C2, . . .) be the sequence of numbers where c\ is the number of initial south 
s steps in Bounce('P), C2 is the number of contiguous west steps that follow the initial run of south 
steps in Bounce(T'), and so forth. 

Example 1.2. The bounce path is indicated by a thick directed line in the following diagrams: 
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For Vi we have Bounce('Pi) = (s, s, w, w, w, s, w, s, s) and bounce('Pi) = (2,3,1,1,2). For V2 we 
have Bounce("P2) = (s, w, w, w, s, s, w, w, w) and bounce("P2) = (1,3,2,3). 

In the appendix we give two characterisations of parallelogram polyominoes. Characterisa- 
tion A.l is in terms of the bounce path of a polyomino. Characterisation |A.5 is in terms of 
the two integer sequences from which a polyomino may be constructed. 

2. The sandpile model on D m , n and recurrent configurations 

In this section we show how recurrent configurations of the sandpile model [9] on K m ^ n which 
has a designated vertex that acts as a sink can be classified in terms of parallelogram polyominoes. 
In what follows we will refer to the the directed bipartite graph K m ^ n that has a designated vertex, 
which we will call the sink, as D m ^ n . The purpose in doing so is to avoid confusion about two 
non-equivalent choices for the sink in K m ^ n (in general if we have many sinks then we merge them 
into a single one). 

For a general (stable) configuration u on F) m n we define a collection of cells fm^niy^) — [0,772] x 
[0, n] C R 2 . A canonical toppling process for checking recurrent configurations is defined (using 
Lemma 2.2) for the graph D m ^ n , and it is then shown that a configuration being recurrent is 
equivalent to f m ,n(u) being a parallelogram polyomino in Para mjn . During this, we show the 
canonical toppling process of a recurrent configuration to be intimately linked with the bounce 
path of the corresponding polyomino f m ,n(u). 

Let -D m ,n be the directed graph D m ^ n (y, E) with vertex set V = {vo, . . . , v m+n -i} and a directed 
edge set 

E = {(vi, Vj) : min(«, j) G {0, . . . , m — 1} and max(i, j) G {m, . . . , m + n — 1}}. 

This is illustrated in Figure [3j We call {vo, . . . , f m _i} and {v m , . . . , v m+n -i} the top vertices and 
bottom vertices of D m>n , respectively. 
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Figure 2 . How the different types of recurrent configurations of Section 3 relate to 
one another. The intersection of the line-filled ellipse with the line-filled rectangular 
region is the collection of "minanz" recurrent configurations. When m = n this is the 
set SqRec n defined and discussed in Section 3. When m = n, the elliptical region 
contained within the minanz region is the collection of square top-heavy minanz 
recurrent configurations SqRec^, also defined and discussed in Section 3. 




Figure 3. The graph D m>n . 

A configuration, or state, x on D m ^ n is an assignment 

x : {1, . . . ,m + n - 1} ->■ N = N U {0}, 

the collection of non-negative integers. The value x\ = x{%) represents the number of grains of sand 
on vertex v Vertex vq is treated as the sink and the number of grains of sand on this vertex is 
generally ignored. There is, however, one benefit to considering the number of grains at the sink 



and this will be mentioned in the proof of Theorem 2.6 

Let eij = 1 if (vi,Vj) S E, and otherwise. Let di be the total out-degree of vertex «j. The 
toppling aspect of the sandpile model is as follows: Suppose x = {x\, . . . ,x m+n -i) is a state on 
D m>n . If Xi > di for some i > 1, then let y = (yi, . . . , y m+n _i) be the state with yi = X{ — di and set 

Vj = x j + e ij f° r an j 7^ This is denoted x A y and means that in configuration x the unstable 
vertex Vi topples leading to the configuration y. 
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One continues this toppling procedure until X{ < di for all % G {1, . . . , m + n — 1}. (The order 
of the topplings does not matter.) Such a state (xi, . . . , x m+n _i) is called stable. Let Stable(D mjn ) 
be the set of stable states on D m ^ n . For example, it is easy to see that |Stable(Z? mri )| — n m ^m n . 

Let us call a configuration u G Stable(D mj?1 ) increasing if u\ < ■ ■ ■ < u m -\ and u m < ■ ■ ■ < 
Um+n-i. Every configuration u G Stable(D mjn ) may be written uniquely as a pair decomp mn (u) = 
(inc mjn (n), sperm m n (n)), where '\nc m ^ n {u) = u' is an increasing configuration, sperm ?nn (n) = ir 
is the lexicographically smallest permutation such that n : {l,...,m — 1} — > {1, . . . ,m — 1}, 
7r : {m, . . . , m + n — 1} — > {m, . . . ,m + n — 1}, and u» = u' n ^ for all i G {1, . . . , m + n — 1}. 

Example 2.1. Consider u = (5,2,1,2,4,2,1,3,2,2,1) G Stable(Z? 5j7 ). VFe /icwe inc 5i7 (u) = 
(1, 2, 2, 5, 1, 1, 2, 2, 2, 3, 4) and sperm 57 (u) = (4, 2, 1, 3, 11, 7, 5, 10, 8, 9, 6).' 

Let o~(x) be the stable state that results from initial state x. A stable state is recurrent if it is 
<t(x) for some x with Xi > di for 1 < i < m + n — 1. Let Rec(D mjn ) be the set of recurrent states on 
D m ^ n and let Rec* \D mtn ) be the set of increasing recurrent states on D mtn . Recurrent states may 
be classified for general directed graphs by the following lemma (see [19] or |9] ) . 

Lemma 2.2. Let D be a directed graph on the vertices Vq, V\, . . . , V n with eij arcs from Vi to Vj, 
and en = 0. Let D be balanced (the in-degree of each vertex equals its out-degree). Then the state 
x = (xi, . . . , x n ) is recurrent iff x = a(x + a') where a' = (eoi, • • • , eo n ). Furthermore, in toppling 
from x + a' to x, every vertex V\, . . . , V n topples exactly once. 

The above lemma applied to the graph D m ^ n gives us that a stable configuration u G Stable(Z) mjn ) 
is recurrent iff the Dhar criterion [9j holds: 

a(ui, . . . ,u m -i, 1 + 1 + ) = («l,... )• (1) 

The order of topplings is unimportant when checking that the Dhar criterion |l]) holds. However, 
for our purposes it will prove useful to fix a canonical toppling process as follows: Let Qi be the 
set of all unstable vertices in the bottom row of D m ^ n (as a result of adding 1 to the height of each 
vertex in this bottom row). Topple all vertices in Q\ and then let P\ be the set of all unstable 
vertices in the top row of D m ^ n . Topple all vertices in Pi and then let Q2 be the set of all unstable 
vertices in the bottom row of D m ^ n . Continue this process until all vertices of D m ^ n are stable. Let 

CanonTop(u) = (Q\, Pi, Q 2 , P2, ■ ■ ■), and 
canontop(u) = (|Qi|, \Pi\, \Q 2 \, \P 2 \, . • •)• 
We note of course that it is possible for CanonTop(u) = canontop(u) = (), the empty sequence. 

Example 2.3. (i) Consider u = (2,0,2,1,0,2) G Stab^L^). When one grain is added to 
every vertex in the bottom row, we have the unstable configuration w = (2,0,3,2,1,3). 
Alternate between toppling all unstable vertices in the bottom row and top row: 

(2, 0, 3, 2, 1, 3) -»■ (4, 2, 0, 2, 1, 0) -> (0, 2, 1, 3, 2, 1) -> (1, 3, 1, 0, 2, 1) G Stable(L> 3 ,4). 

Since a(w) = (1, 3, 1,0, 2, 1) 7^ u, u is not recurrent. 

Vertices v% and v% were first to topple in the bottom row, thus Q\ = {^3,^6}- These 
caused vertex v\ in the top row to become unstable, so P\ = {v\}. Once vertex v\ toppled, 
it caused vertex in the bottom row to become unstable, so Q 2 = {v^}. When vertex 
topples it does not cause any vertices in the top row to become unstable and the toppling 
process now ends. We have CanonTop(n) = ({3, 6}, {1}, {4}). 
(ii) Consider u = (0, 2, 1, 2, 1, 2) G Stable(L>3 5 4). Adding one grain to every vertex in the bottom 
row gives w = (0,2,2,3,2,3), which topples as follows: (0,2,2,3,2,3) — > (2,4,2,0,2,0) — > 
(2,0,3,1,3,1) (4,2,0,1,0,1) (0,2,1,2,1,2). Since a(w) = u we have that u G 
Rec(Z?3 5 4). The toppling process is CanonTop(n) = ({4, 6}, {2}, {3, 5}, {1}). 
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In the next theorem we will classify the set Rec(-D mjn ), however some terminology is required 
first. 

Definition 2.4. Given u = (m, , u m +n~i) G Stable(L> mjn ) with CanonTop(-u) = {Qi, Pi,Q2, P2, • • 

and inc m n (u) = (ai, . . . , a m _i, 61, ... , b n ), define f m n to be the collection of cells in the plane given 
by; 

(m— 1 \ / n \ 

[m-l,m] x [0,n]U |J [t-l,i] x [0,1 + a.] J n f (J[0,l + 6i] x [i-l,*]J . 

Let us note that f m ,n{u) is the intersection of the Young diagram whose corner at (m,0) has 
heights (n, 1 + a m _i, . . . , 1 + a±) from right to left, and the Young diagram whose corner at (0, n) 
has widths (1 + b n , . . . , 1 + 61) from top to bottom. 

Example 2.5. Given u = (2, 1, 2, 0, 0, 2, 6, 1, 5, 1) G Stable(Z>7 ) 4), we /icwe 

inc 7i4 (V) = (01, . . . , a 6 , 61, . . . , 64) = (0, 0, 1, 2, 2, 2, 1, 1, 5, 6). 
The collections of cells f7,4(u) is 
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Since fniu) Para 74, the next theorem proves that u is not recurrent. 
Theorem 2.6. Let u = (u±, . . . ,u m+n _i) G Stable(-D mjn ). 

(a) u G Rec(An,n) iff f m ,n{ u ) G Para 

m,n • 

(b) If u £ Rec(D rrhn ) then bounce(/ mjn (u)) = canontop(u). 



Proof. Let u = (ui, . . . ,u m+ri _i) G Stable(D mjra ). It is a consequence of Lemma 2.2 that in testing 
whether u is a recurrent configuration, every vertex in {^i, . . • 5 f^-fn— l 

} topples exactly once during 

the process. The configuration (u\, . . . , u m -i, 1 + u m , . . . , 1 + Um+n-i) may be interpreted as the 
result of the toppling of the sink vq in the configuration u even if the height of this vertex is not 
specified. If every other vertex in {v\, . . . , v m+n -\} topples exactly once, then the height of vertex 

v% will be Ui + (^2j — n = Ui for 1 < i < m — 1, and 1 + m + [Y^j^Q ej,^ — m = Ui for all 
m<i<m + n — 1. The order of topplings is irrelevant, and in what follows we will adhere to the 
canonical toppling process mentioned earlier in this section. 

The state u is recurrent iff u' = inc m>n («) is recurrent, since it is simply a reordering of heights 
within the top and bottom rows of the graph. If CanonTop(u) = (A±, A2, ■ ■ •) and CanonTop(u') = 
(A^, A' 2 , . . .) then Ai = {vr(j) : j G A' { } for all i > 1, where ir = sperm m n (n). Consequently, we 
have canontop m n (» = canontop m n (-u')- 

Dhar's criterion Q applied to v! is equivalent to there existing non-empty ordered partitions 

P = {Pi...,P' k ) andg = (g' 1 ,...,Q' fc+1 ) 
of the sets {vi, . . . , v m -i} and {v m , . . . , v m+n -\}, respectively, with the following properties: 
(i') 1 + u'i = m for all Vi G Q' v 

(ii') u[ + q + . . .qj > n > u ■ + q + . . . + qj_i for all V{ G Fj, 1 < j < k. 
(hi') 1 + u'i +po + . . . +pj > m > 1 + u'i +po + . . . + pj-i for all Vi G Qj + i, 1 < j < k. 

where p» = |P/|, % = and p = 5o = 0. 

The understanding above is that Q[ is the set of vertices on the bottom row that topple first, 
i.e. Q\ = {vi : 1 + u[ > m, 1 < i < m — 1}. Then we topple the vertices in the top row that have 
become unstable due to grains being added from the previous topplings in the bottom row, i.e. 
P[ = {vi : 1 + u' i + q\ > n, m<i<m + n — 1}. Next we topple those vertices in the bottom row 
that were stable before the topplings of the vertices P[, but have consequently become unstable, 
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i.e. Q' 2 = {vi : 1 + U{ < n, 1 + v! i + px > n, m < i < m + n — 1}, and so on. The block topplings 
happen via Q[ — > P[ — > Q' 2 — > P 2 — > Q3 • • • — > Q'k+v All °f these sets except for Q' k+1 must be 
non-empty. 

Since the initial heights of the vertices in the top row of u' are weakly increasing, the vertices 
may be partitioned into contiguous blocks according to the 'time' at which they topple. The same 
is true of the bottom row. Indeed po, . . . ,p k and qo,---iQk are such that 

Vi £ P'j % G {m - (p + ■ ■ ■ + pj), ■ ■ ■ , m - 1 - (p + . . . + Pj-i)} 

Vi G Qj 44> i G {m + n - (q + . . . + qj), ...,m + n — 1 — (q + ■ ■ ■ + qj-i)}- 

with 

po + . . . + p k = m-l 
qo + • • • + qk+i = n 
P!,...,p k ,qi,...,q k > 0. 

Note that it is possible that qt+i = 0. Let us now change variables by setting a% = for all 1 < i < 
m-l,bi = u' m _ 1+i for all 1 < i < n, Xi = m - 1 - (p + ■ • • +Pk-i) and m = n - (q + . . . + qk+i-i) 
for all suitably defined i. Conditions (i')-(iii') are equivalent to the following: 

There exist numbers = xq < x\ < ■ ■ ■ < Xk = m — 1 and = yo < y\ < ■ ■ ■ < yk+i = n such 
that: 

(i") 1 + bt = m for all 1 + y k < i < y k +\ = n. 

(ii") 1 + Oj G {1 + yj+i, • • • , yj+2} for all 1 + Xj < i < Xj + i and < j < k — 1. 
(iii") 1 + 6j G {1 + Xj, . . . , Xj+i} for all 1 + yj < i < yj+i] and < j < k — 1. 

By defining x^+i = m, conditions (i") and (iii") above may be more compactly written as 

1 + bi G {1 + Xj, . . . , Xj+i} for all 1 + yj < i < yj+i and < j < k. 

We have now shown that given u G Stable(D m . n ), the configuration u G Rec(Z? mjn ) iff there exist 
sequences (x ,..., x k+1 ) and (yo,..., y k +i) with 

(i) = xo < x\ < • • • < x k = m — 1 < x k+ \ = m, 

(ii) = y < y x < y 2 < • • • < y k +i = n, 

(iii) 1 + a,i G {1 + yj+x, ■ ■ ■ , yj+2} for all 1 + Xj < i < Xj + x and < j < k — 1, 

(iv) 1 + h G {1 + Xj, . . . , Xj + x} for all 1 + yj < i < yj+x and < j < k. 



From Characterization A.l these conditions hold precisely when f m ,n(v) G Para mifl . This completes 
part (a). 

The bounce path of such a parallelogram polyomino is (yk+x — Vk,x k — x k -x,y k — Vk-i,Xk-i — 
x k _2, ■ ■ •) = (qi,Pi,q2,P2, ■ ■ •)• Since the vertices topple in the order Qx —> Px —> Q k +x (note 

that Q k +x can be empty) we have that canontop(u) = (qx,Pi,q2,P2, ■ ■ •) = bounce(/ min (n)), which 
concludes part (b). □ 

One may give a bijection from Rec(D myn ) to the set of all parallelogram polyominoes whose 
bounce paths have been decorated with an ordered set partition as follows: 

Associate to every parallelogram polyomino a pair of ordered set partitions whose number of parts 
and size are related to a path within the polyomino that traverses the polyomino by 'bouncing off' 
edges and corners. 

If bounce('P) = (cx,C2, ■ ■ •), then let us write 

ebounce^) = (02, C4, . . .) and obounce^) = (cj, C3, . . .). 

Given a set A we will write X = (X\,X2, ■ ■ ■) 1= A to mean that X is an ordered set partition 
of A. We will write type(X) = (|-Xi|, \X 2 \, . . .), the sequence of sizes of the ordered sets in the 
partition. For example, if X = ({3, 4, 2}, {1, 5}, {6, 7}) N {1, . . . , 7}, then type(X) = (3, 2, 2). 
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Let SetPara min be the set of all triples (V,A,B) where V € Para min , A = (Ai, A2, ■ ■ .) is an 
ordered set partition of {1, . . . ,m — 1} where Ai is associated with the ith run of west steps in 
Bounce('P), and B = {B\, B2, ■ ■ ■) is an ordered set partition of {m, . . . ,m + n — 1} where Bi is 
associated with the ith run of south steps in Bounce('P). More formally, 



SetPara r 



U U 



u 



(V,A,B). 



VePara m ,n A(={l,...,m-l} B(={m,...,m+n-l} 

type( A)— ebounce('P) type(_B)=obounce{'P) 



Note that if (V,A = (A 1 ,A 2 ,...),B = (B 1 ,B 2 ,...)) G SetPara m , n 
then In = \AA and I21-1 = \P%\ f° r au £ > 1. 



Example 2.7. Let V be V2 in Example 1.2, Let A 
Then (V,A,B) € SetPara7 ; 3 and may be represented in the following way: 



and bounce(P) 

({1,4, 5}, {2, 3, 6}) and B = 



= (h,l2,-..), 

({8}, {7, 9}). 



{2, 3, 6} {1,4, 5} 



(V,A,B) 



























<— 








<- 






> 












<■■ 

r 











{8} 
{7,9} 



Given u £ Stable(Z) min ) with CanonTop(u) = (Qi, Pi, Q2, P2, ■ ■ ■), let 

n(u) = (/ m , n (u), (Pi,P 2 , . . .), (Ql, Q2, • • •))• 



This construction, combined with parts (a) and (b) of Theorem 2.6, gives: 

Corollary 2.8. Q, : Rec(D mjn ) — > SetPara mj „ is a bisection. 

It is well-known from the literature on the sandpile model that recurrent configurations on a 
graph G are in one-to-one correspondence with spanning trees of G. Thus the number of recurrent 
configurations on D m , n is the number of spanning trees of K m ^ n , which is m 



n— l„m— 1 



Corollary 2.9. For allm,n > 1, |Rec(D 
where Nara(a, b) — 1 ■" 



m n 1 ri m 1 



and \Rec* (D m ^ n )\ = Nara(m+n— 1, m), 
aKbJKb-i) are the Narayana numbers |22| . 



3. Special types of minimal recurrent configurations 

In this section we examine several different types of minimal recurrent configurations in Rec(-D m)n ) 
under the map / that was given in the previous section. As was mentioned before, the relation- 
ship between these different collections is summarised in Figure [2| Recurrent configurations which 
contain no empty vertices are a natural class to examine from a physical viewpoint. However, if a 
recurrent configuration is minimal, then there must be at least one vertex that is empty (this will 
be shown in Lemma 3.4). So in what follows we will be looking at minimal recurrent configura- 



tions that are almost non-empty, which we will call minanz configurations, in the sense that one 
designated vertex on the bottom row is allowed to be empty. 

In Section 3.2, we look at configurations that have the same number of vertices on the top 
row as on the bottom row. We call such configurations square configurations, i.e. those for which 
m = n. Square configurations in Rec(D nin ) that are minanz are shown to correspond via a mapping 
to a matrix analogue of set partitions. After this we consider a subset of the aforementioned 
square configurations, which we will call top-heavy. The image of these configurations under (f> are 



percisely those matrices which are upper-triangular (Theorem 3.7). This leads to a connection with 
(2+2)-free partially ordered sets, which is explained towards the end of this section (Theorem 3.12 ). 
Given u € Rec(D m ^ n ), let 



level(u) 



ui + .. .+u r , 



mm 



I). 
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Define the equivalence relation ~ as follows: for u,u' G Rec(Z) mjra ), we write u ~ v! if inc mjn (n) = 
inc m , j?1 (u / ). The quotient set Rec(.D m]ri )/ ~ is in one-to-one correspondence with Rec* (D mn ), the 
set of all increasing configurations of Rec{D m<n ). 

Theorem 3.1. If u £ Rec(D min ), then level(u) = area(/ min (u)) — (m + n — 1). 
Proof. Let u G Rec(£> m n ). Then we have that 



m—1 m—1 

n 

i=l i=l 



+ ^2 (1 + u%) = n + ^2 ( l + a ») = area (/m,»M) + larea(/ mi „(n)) 



and 



m+n— 1 



^2 (1 + Uj) = ^(1 + bj) = area(/ m>n (it)) + uarea(/ mj „(u)). 

j=m j=l 

Adding both of these gives 

m+n— 1 

n+ ^2 ( 1 + u i) = uarea(/ m , n (n)) + 2area(/ m , n (V)) + larea(/ TO)n (u)) 

i=l 

= mn + area(/ min (u)), 

which gives level(u) = area(/ mi „(n)) — (m + n — 1). □ 

3.1. Minimal and minanz recurrent configurations. A configuration ti G Rec(Z? m n ) for which 
level (u) is as small as possible, i.e. 0, is called minimal. From Theorem |3.1[ level (u) = iff 
area(/ mi „(n)) = m + n — 1, and this happens iff f m ,n(u) G Ribbon mjn . 

Theorem 3.2. Let u G Rec(D mtn ) . The following are equivalent: (i) u is minimal, (ii) f m ,n(u) G 
Ribbon mjn , and (Hi) level(u) = 0. 

The set of all different minimal configurations {u G Rec*(D mi „) : u is minimal} is in one-to- 
one correspondence with the set of all ribbon polyominoes Ribbon m n . Every ribbon polyomino is 
uniquely characterised by a lattice path that goes through all the centres of its cells. Such a path 
goes from (1/2, 1/2) to (m — 1/2, n — 1/2) and takes unit north and east steps. 

Corollary 3.3. \{u G Rec*(D mjn ) : u is minimal}\ = ( m ^ 2 ) ■ 

We call a configuration u G Rec(D m)n ) almost non-zero if 

U\ , . • • , Um— 1 j ^m+1 > • • • > ^m+n— 1 > 0. 

Configurations that are both minimal and almost non-zero will be called minanz configurations. 

Lemma 3.4. // u G Rec(Z) mjn ) is minanz, then 

(i) u m = ; 

(ii) the cells [0,1] x [0,1], [0,1] x [1,2] and [1,2] x [1,2] are in f m ,n(u), and 

(iii) Bounce(/ mjn (u)) has suffix ws. 

Proof. If u G Rec(D m ^ n ) is minimal, then V = fm,n{u) is a ribbon polyomino. A ribbon polyomino 
contains either the cell [1, 2] x [0, 1] or the cell [0, 1] x [1, 2], but not both. Since U\, . . . , u m _i > 0, 
we have that 1 + u±, . . . , 1 + u m -\ > 2. Thus [0, 1] x [1, 2] is in V and [1, 2] x [0, 1] is not, which 
implies the smallest entry of u m , . . . , u m+n -\ is 0. Since u m+ i, . . . , ii m + n _i > 0, we must have 
u m = and also that [1, 2] x [1, 2] G V. 

For a general u G Rec(Z) mjn ), Bounce(/ mi „(u)) may end with either of the four suffixes: ss, 
sw, ws, or ww. Suppose that Bounce(/ m>n (n)) ends in ww. One must then have u m > 0, which 
contradicts (i) above. If Bounce(/ m]?1 (u)) ends with sw, then consider the third last step. The suffix 
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of Bounce(/ mjn (u)) is either (a) ssw or (b) wsw. Neither (a) nor (b) can be bounce paths in the 
final cells of Bounce(f m ^ n (u)) given in (ii). Finally, suppose that Bounce(/ min (ii)) ends in ss. The 
two steps preceding ss must then be ww. Considering the position of the leftmost cell at height 
[2, 3] shows that no such path exists with suffix wwss. □ 

For a general graph G, minimal recurrent configurations on G are in bijection with acyclic 
orientations of G whose unique sink is a distinguished vertex, as is the sink in the sandpile model 
we are considering. Our restriction to minanz configurations corresponds to acyclic orientations 
where in addition there is a unique fixed source that shares an edge with the sink. (See Gioan and 
Las Vergnas [16"].) 

It is a consequence of (iii) and the identification of the bounce path with canontop(ii) that if 
u G Rec(D m)n ) is minanz with CanonTop(u) = (Qi, Pi, ... , Qk+i), then Qk+i = {m}. 

Corollary 3.5. \{u G Rec*(.D TO) „) : u is minanz}\ = ( m ^™2 4 ) • 
Next we consider the special case for when m = n. 

3.2. Square minanz and square top-heavy minanz recurrent configurations. Let BiComp n 

be the set of square matrices whose entries partition the set {1, . . . , n} and having no rows or 
columns consisting of only empty sets. Let BiComp,^ be the collection of upper triangular matrices 
in BiComp n . Define 

SqRec n = {«£ Rec{D n ^ n ) : u is minanz}. 

Suppose that u G SqRec n with CanonTop(n) = {Qi, Pi, . . . , Qk+i = { n })- If vertex v x is such 
that x G Pi U Qi then we say that v x is in the ith wave of n and denote this by wave u (v x ) = i. Let 
us define 

SqRec^ = {u G SqRec n : \Nave u (v n + x ) < wave u (t> x ) for all 1 < x < n}. 

In words, SqRec^ is the collection of all square minanz configurations whose canonical topplings 
have the property that bottom vertices topple on or before the wave which topples the vertex 
directly above. 

From the definition of SqRec^, and using the fact that /(SqRec^) C Ribbon„ in , the image 
/(SqRec^) is precisely the class of ribbon polyominoes whose cells/squares are such that none of 
their centres lie beneath the line y = x in the plane. 

Given u G Rec(D n ^ n ) is minanz with CanonTop(-u) = (Qi, Pi, . ■ ■ , Qk+i = { n }), define Q\ = Qi—n 
for all i. So (Q[, . . . , Q' k ) and (Pi, . . . , Pk) are both ordered set partitions of {1, . . . , n — 1}. Let 
4>{u) be the k x k matrix M with Mij = Pi n Q'j- Equivalently, we have that x G My if x G Pi and 
n + x G Qj for all x G {1, . . . , n — 1}. 

Example 3.6. (i) Consider u = (4, 5, 6, 1, 4, 5, 4, 0, 7, 1, 1, 4, 6, 1, 7) G SqRec 8 . We have 
CanonTop(n) = ({9, 15}, {3}, {13}, {2, 6}, {12}, {1, 5, 7}, {10, 11, 14}, {4}, {8}) 
so that 

(Pl,P 2 ,P 3 ,P A ) = ({3}, {2, 6}, {1,5, 7}, {4}) 
(Qi,Q2,Q 3 ,Q±) = ({9, 15}, {13}, {12}, {10, 11, 14}) 
=> (Qi, C&, Qa, Qi) = ({1, 7}, {5}, {4}, {2, 3, 6}). 



Thus 



( {3} \ 

{2,6} 

{1,7} {5} 

\ {4} J 
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(ii) Consider u = (4, 3, 4, 1,0, 2,1,4, 1) G SqRec 5 . We have 

CanonTop(u) = ({8}, {1, 3}, {6}, {2}, {7, 9}, {4}, {5}). 

From this (P U P 2 , P 3 ) = ({1, 3}, {2}, {4}), (Q 1; Q 2 , Q 3 ) = ({8}, {6}, {7, 9}) and so (Q[,Q' 2 , Q' 3 ) = 
({3},{1},{2,4}). Thus 

/{3} {1} \ 
<j)(u) = {2} . 
V ^ {4} ) 

Theorem 3.7. (i) <fi : SqRec n — > BiComp n is a bijection. (ii) (f> : SqRec^ — > BiComp n Tq is a 

bijection. 

Proof. For (i) let u G SqRec n . We have CanonTop(u) = (Qi, -Pi, . . • , Qk+i = {n})- From the nature 
of the toppling process (Q' 1; . . . , Q' k ) (= {1, . . . , n— 1} and (Pi, . . . , P&) N {1, . . . , n— 1}. Since every 
element x G {1, . . . , n — 1} appears in exactly one of the P sets, Pi say, and appears in exactly one 
of the Q' sets, Q'j say, we have that x G Mij = Pi n Q'-. There can be no row or column of empty 
sets for this reason. Thus M = <f>(u) G BiComp n . 

Let us now define a function [i on BiComp n . Suppose that M G BiComp n with dim(M) = k. For 
all ie{l,...,k}, let 

Pi = Ma U . . . U Af ifc 
Q'i = M u U . . . U M w . 

Define = Q'j + n for all i G {1, . . . , k] and Qk+i = {n}. Let Po = qo = 0, ^ = |Pj| and % = |Qj|. 
Let u = /i(M) be the configuration where u n = 0, 

u n+x = n — 1 — (p + • • • + Pi-i) for all x e Qi and i G {1, . . . , &} 
u x = n — 1 — (go + • • • + Qi) for all x G Pi and i G {1, . . . , k}. 

The value level (u) is minimal. Since Qfc+i = {n}, u n = n — 1 — (po + ■ ■ ■ + Pk) = 0, and Ui > for 
all other i, hence u is almost non-zero. Thus u = /x(M) G SqRec ra . It is straightforward to check 
that (j){iJ.{M)) = M and /j,(4>(u)) = u for all M G BiComp n and u G SqRec n , respectively. Thus 
/i = 4>~ l and : SqRec n — > BiComp n is a bijection. 

Part (ii): from part (i) we have that (/>(SqRec^) C BiComp n . Let u G SqRec n , M = (p(u) 
and suppose that dim(M) = k. The matrix M is upper-traingular (i.e. M G BiComp„) iff Q\ C 
Pi U • • • U Pi for all 1 < i < fc. This condition is equivalent to 'wave u (u n+:r ) = i ^> wave u (v x ) < i 
for all 1 < i < and 1 < x < n — 1', which in turn is equivalent to 1 wave 11 (f n + a; ) < wave u (v x ) for 
all 1 < x < n - 1. □ 

Corollary 3.8. \{u G Rec(P> nj „) : u is minanz}\ = ^2kS(n — l,k) , where S(n,k) are the 
Stirling numbers of the 2nd kind. 

As we mentioned in the introduction, the matrices BiComp n ^ were the subject of the two papers [51 
[TO] which related these matrices to (2+2)-free partially ordered sets (posets) on the set {1, . . . ,n} 
(also known as interval orders). There are several other rich connections between these objects and 
labelled Stoimenow matchings, ascent sequences and pattern avoiding permutations. The remainder 
of this section is dedicated to showing a direct link between (2+2)-free posets and the recurrent 
configurations of SqRec^. This link was the original motivation behind this paper. 

Let P = (P, -<p) be a partially ordered set on {1, . . . , n}. We denote by P* the dual poset of P. 
Given x G P, let D(x) = {y G P : y -< x} be the down set of x. A defining property of (2+2)-free 
posets is the following: a poset P is (2+2)-free iff the set of all down-sets of P may be linearly ordered 
by inclusion. Let Posets n be the set of all (2+2)-free posets on the set {1, . . . , n}. Every P G Posets„ 
can be uniquely represented as a sequence of sets P = (Dq(P), . . . , Dk~\{P) \ Lq(P), . . . , Lk-i(P)) 
where = P>o(P) C -Di(P) C • • • C £> fc _i(P) and (L (P), ... , L k -\{P)) is an ordered set partition 
of {1, . . . , n} such that if x G Lj(P) then D(x) = A(P)- 
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Definition 3.9. Given M G BiComp^, let P = (P, <p) = mtp(M) be the partially ordered set on 
{1, . . . ,n} where x <p y iff there exist i,j, i',f with x G My, y G Mjy and j < i! . 

We refer the reader to the papers [51 [10] if there is any confusion surrounding the definitions or 
terminology. 

/{3} {7,2} \ 

Example 3.10. Let M = \ {5} {1} . Then P = (P, ^) is tfte pose£ on {1,2,..., 7} 

\ {4,6}/ 
with the following relations: 3 ~<p 1, 4, 5, 6 and 2,5,7 -<p 4, 6. 

4 6 




Theorem 3.11 (Dukes et al. [10]). Let Posets n be the set of all (2+2)-free posets on {1, . . . , n}. 
Then mtp : BiComp^ — > Posets„ is a bijection. 

In the next theorem we see that the composition of the functions <p and mtp is such that the 
recurrent configuration can be easily read from the level structure of the associated poset. 

Theorem 3.12. Let u G SqRec^ with CanonTop(n) = (Qi, Pi, . . . , Q k +i = {n}) and (P, ■<) = 
mtp((/)(n)). Then 

u x+n = \D j+1 (P*)\ for all x G Lj(P*) and j G {0, . . . ,k - 1} 
u x = n — \Dj + i(P)\ for all x G Lj(P) and j G {0, . . . , k — 1}. 

Proof. Let u G SqRec^ with CanonTop(n) = (Qi, Pi, ■ ■ ■ ,Qk+i = { n })- Let M = cf)(u) and 
(P, -<) = mtp(M). We have dim(M) = k, Pi is the union of entries in row i of M, and Qi is the 
union of entries in column i of M with all entries increased by n, for all i G {1, . . . , k}. The poset 
P = (P, <p) is specified by; 

P = (A,(P), • • • , P fc -i(P); MP), • • • > L k -x(P)). 

The dual P* of P is given by 

P* = (A)(P*), ■ ■ ■ , D k ^(P*); L (P*), . . . , L fe -i(P*)) 

where L fe _i_i(P*) = A+i(P)\A(P) and D^^P*) = L i+ i(P) U ... U L A -i(P) for all i G 
{0, . . . , k — 1}. We use the convention Pfc(P) = {1, . . . , n — 1}. 



From Theorem 3.7 we have u n = 0, 

u x+n = n — 1 — (po + . . . + Pi-i) for all x £ Qi and i € {1, . . . , /c}, and 
u x = n — 1 — (qo + ■ ■ ■ + qi) for all x G Pj and i G {1, . . . , /c}. 

For all i 6 {1, ... , k}, Pi = L^P) and Qi = A(P)\A-i(P) = L k ^(P*). For alH G {1, . . . , fc}, 

n - 1 - (p + . . . + Pi-i) = pi + ...+pk 

= |L i _ 1 (P)U...UL fc _ 1 (P)| 

= |£> fc _i_<(P*)| 

and 

n - 1 - (g + • • • + qi) = Qi+i + ... + % + Qk+i 

= \Q i+1 U...UQ k \ + q k+1 

= |(A+i(P)\A(P)) u . . . U ( J D fe (P)\D fc _ 1 (P))| + 1 

= \D k (P)\Di(P)\ + 1 

= n-l-|A(P)| + l = n-\Di(P)\. 
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Example 3.13. Consider the poset P given in Example 3.10 , The level- and down-sets of this 
poset are illustrated in the following diagram: 



P = 








D 3 (P) 


= {1,2,3,4,5,6,7} 


L 2 (P) 


= {4,6} 


D 2 (P) 


= {2,3,5,7} 


MP) 


= {1,5} 


Di(P) 


= {3} 


MP) 


= {2,3,7} 


Do(P) 


= 



From Theorem 3.12 we have that u x = \Dj + i(P)\ for all x G Lj(P) and j G {0,1,2}, hence 
U2 = u 3 = u 7 = 8 - 1 = 7, ui = u 5 = 8 - 4 = 4, u 3 = 7 => ...,u 7 ) = (4,7,7, 1,4, 1,7). 
The dual P* of P, together with its level- and down-sets: 

D 3 (P*) = {1,2,3,4,5,6,7} 
£ 2 (P*) = {3} D 2 (P*) ={1,4,5,6} 

- . r, ..... f -2 .7 Li(P*) = {2,5,7} D 1 (P*) = {4, 6} 

L (P*) = {1,4,6} J D o ( J P*)=0 




4 6 



From Theorem 3.12 we have thatu x = \Dj + \{P*)\ for allx £ Lj(P*) and j £ {0, 1, 2}, hence us+i = 

u 8+7 = 4 and u 8 +3 = 7. So (u 9 , . . . , u 15 ) = (2, 4, 7, 2, 4, 2,4). TTie 



-u-8+4 = w 8+6 = 2, u 8+2 = ^8+5 
configuration is u = (4, 7, 7, 1,4, 1, 7,0,2,4,7,2,4,2,4). 

We end this section with the following conjecture concerning square non-zero recurrent configu- 
rations. Is there a combinatorial explanation for the relationship to walks in the plane? (See |2,3l 
Sequence A145600].) 

Conjecture 3.14. Let a n = {u £ Rec*(£> n , n ) : u u . . . , n m+n _i > 0}| . Then a n = ^ ( 2n ~ 2 ) ( n 2 ™ 2 ) ; 
which is the number of walks from (0,0) to (0, 1) that remain in the upper half-plane (y > 0) using 
2n — 3 units steps {n, s, e, w} 



4. q, i-NARAYANA POLYNOMIALS AND THEIR SYMMETRY 

In this section we will introduce a polynomial that we call the g, t-Narayana polynomial. The 
polynomial is the generating function for the bistatistic (area, parabounce) on the set of parallelo- 
gram polyominoes. In terms of recurrent configurations of the sandpile model, the area statistic 
is proportional to the sum of the heights of the piles, and the parabounce statistic is related to 
the canonical toppling process. When viewed in the parallelogram polyomino, the bounce path is 
almost identical to the bounce path that Haglund defined for Dyck paths [18J. Our polynomial is a 
natural extension of the area and bounce path statistics to the class of parallelogram polyominoes. 
More will be discussed about this in the subsections that follow. 

For any polyomino V G Para m>n , define its bounce weight to be 

k 

parabounce('P) = "}Ji/2\cj, 
i=i 



where bounce('P) = [c\, . . . c^). 
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This weight may also be described by summing weights on each step of the bounce path. The 
initial step has a weight of 1, and the weight of subsequent steps is incremented by 1 after each 
turn on the upper boundary path of the polyomino. Consider the following two polyominoes: 



Vx 













4- 












& 


2 
t 
















3 












f 


3 




3 

r 










4 
t 
















5 


















P 















e <Para 9 7 ; 



v 2 









,1 


1 


1 




4- 


2 


2p 











G <Para s ,2- 



The bounce weights are 

parabounce(T'i) = (1) + (1 + 1 + 1) + (2) + (2) + (3 + 3) + (3 + 3 + 3) + (4) + (4) + (5 + 5) 

parabounce(7> 2 ) = (1) + (1 + 1 + 1 + 1) + (2) + (2 + 2 + 2) = 13. 

The distribution of the bistatistic (area, parabounce) on polyominoes in Para mi „ is represented by 
the generating function 

F m , n (g,t)= g area(P) t parabounce(P) 



41, and 



"PePara™ „ 



We call these polynomials q, t-Narayana polynomials because {F mtn (l, l)} m , n >i are the Narayana 
numbers that were mentioned in Corollary 2.9 The following conjecture has been verified for all 
pairs (m, n) with m < 11 and n G N. 

Conjecture 4.1. For all positive integers m andn, the distribution of the bistatistic (area, parabounce) 
on polyominoes in Para mn is symmetric: 

-^m,n(?)^) — F rnn (t, 0) • 

Further to this we posit another symmetry (which has been checked for all pairs (m, n) with 
max(m, n) < 9). 

Conjecture 4.2. For all positive integers m and n we have F mn (q,t) = F n ^ m {q,t). 

Symmetry along the main diagonal of parallelogram polyominoes provides a bijection from poly- 
ominoes of Para m n to Pa ra n m that preserves the area statistic. This mapping gives the following 



special case of Conjecture 4.2 



Theorem 4.3. For all positive integers m and n we have F min (q, 1) = F nm (q, 1). 



In the next two subsections we will present some arguments that support Conjecture 4.1 In 



the final subsection we will discuss the relationship of Conjecture 4.1 to the sandpile model and a 
similarly symmetric bistatistic on Dyck paths that was introduced by Haglund. 

4.1. Proofs of the conjectures for small values of min(m, n). In this subsection we show 
that Conjectures 4.1 and 4.2 are true for all pairs (m,n) with m = 2. The same method can be 
used to show the conjectures are true for other similarly small values of one of the parameters by 
using regular expressions. These computations will not be detailed in this paper but we list the 
first few polynomials in the appendix, from which symmetry is apparent. Conjecture 4.1 has been 
checked up to min(m, n) = 11 and n £ N using generating functions resulting from the classical 
transfer-matrix method. 

We consider polyominoes in Para2, n for any n > 1. Define the generating function for the 
bistatistic of all these polyominoes, adding a variable z to record the height n of polyominoes: 

F 2 +(q,t;z) = Y / F 2 , n (q,t)z n . 
n>l 
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Theorem 4.4. F 2 , ic (q,t;z) = t 2 q 2 z/{\ - qtz)(l - q 2 tz){\ - qt 2 z). 

Proof. Every polyomino P € Para2,„ is uniquely encoded as a regular expression a*b*cd* by reading 
the polyomino from top to bottom and assigning the different rows the following expressions: 

( a )* ( 3j )* 

( b r ( m r 
( d r («□ r 

In order to count these expressions in terms of the parameters mentioned in the definition of 
F, we weigh the expressions as follows: w(a) = qt is the weight of a row with only the rightmost 
cell occupied, w(b) = q 2 t is the weight of a row with both cells occupied and a vertical step of the 
bounce path between them, w(c) = q 2 t 2 is the weight of the row where the bounce path changes 
from vertical to horizontal, w(d) = qt 2 is the weight of a row with only the leftmost cell and the 
bounce path on a vertical trajectory towards the origin. The straightforward translation of the 
regular expression into a rational generating series gives 

- t2q2z 
F2 ^ t]Z > ~ (l-qtz)(l-q 2 tz)(l-qt 2 zY 

The involution that exchanges the letters b and d, and then exchanges the blocks b* and d* , proves 
the symmetry of i*2,ra(<7; t). □ 

Theorem 4.5. F^ 2 (q,t;z) = t 2 q 2 z/{\ - qtz)(l - q 2 tz){\ - qt 2 z). 

Proof. The case of polyominoes of Para m> 2 is similar to Para2, n , although slightly more complicated 
since it requires a discussion on the height of the first turn in the bounce path. This discussion 
leads to the two cases of the column-by-column decomposition of those polyominoes (details left to 
the reader): 

a* b c* e = g*g*>g and r 9 * h = 

The weights of the expressions as follows: w(a) = qt 2 , w(b) = q 2 t 3 , w(c) = w(g) = q 2 t, w(d) = 
w(e) = w(f) = qt, and w(h) = q 2 t 2 . The translation of the two regular expressions above into 
rational series gives 

F *^ Z) = T^z qH " z -T^zT^z qtz + T^zT^ 

t 2 q 2 z 



(1 — qtz){\ — q 2 tz)(l — qt 2 z) ' 

We may also describe explicitly an involution exchanging the two parameters on the words of this 
regular expression. The equality of weights 

w(be) = w(b)w{e) = q 2 t 3 .qt = q 3 t 4 = qt 2 .q 2 t 2 = w(a)w(h) = w{ah) 

shows that 

w{a*bc*d*e) = w{a + c*d*h), 

with a bijection from a*bc*d*e to a + c*d*h that consists of replacing the letter b by a and the letter 
e by h. In addition, w(g) = q 2 t = w(c) and w(f) = qt = w(d) implies that 

w(f*g*h) = w(c*d*h), 
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with a bijection from f*g*h to c*d*h that consists of replacing the letters / by d and the letters g 
by c, followed by exchanging the order of the blocks of letters d* and c*. Subsequently, one has 

a + c*d*h + c*d*h = a*c*d*h 

and the resulting language is similar to the case of -F2,*> except that now it is the exchange of c and 
a that defines the involution swapping the two statistics. □ 



Corollary 4.6. Conjectures 4-1 an d 4-2 are true for all pairs (m,n) with min(m, n) 
A similar analysis shows that 



] 3 tz 1 — qt 3 2 
See Appendix |A.2| for higher values. 

4.2. Proof of Conjecture |4.1 for minimal values of one of the two statistics. 

Theorem 4.7. In Para m) „, for anyi, there are as many poly ominoes whose bistatistic (area, parabounce) 
is (m + n — as there are polyominoes whose bistatistic (area, parabounce) is (i,m + n — 1). 

In order to prove this theorem, let us assign names to the two sets of polyominoes with which 
we are dealing. The proof will then consist of a bijection T that switches the area and parabounce 
statistics. Let 

A m ,n = {V € Para TO) „ : parabounce(:P) = m + n - 1}, 
B m ,n = {V G Para m>n : area(P) = m + n - 1}. 

Proof. Let 

A k = {[i-l,i] x : i + j - 1 = fc and i,j G N }. 

For any P G Para m)n and all A; € {1, . . . , m + n — 1}, let = dk{V) = \VC\ Afc|, the number of cells 
of "P on its A;th anti-diagonal. Define the vector of all these values diaglen("P) = {d\, . . . d m+n -\) 

Suppose that V £ A m ^ n . We now define an operation T on A m ^ n . Polyominoes in A nitn all have 
the same bounce path that is a sequence of n south steps, followed by a sequence of m — 1 west 
steps. These polyominoes are exactly those which are defined by their upper paths, since their 
bottom paths are necessarily e m n ra . 

The sequence of diagonal lengths of polyominoes in A m ^ n describes this upper path: for i G 
{2, ...,m}, di-i < di (resp. di-\ = di) represents a north (resp. east) step, whereas for i £ 
{m + 1, . . . , m + n}, = di (resp. tfj_i > di) represents a north (resp. east) step. 

The position m is special in the diagonal length sequence since it corresponds to the turn of the 
bottom path: indeed d±, . . .d m is weakly increasing while d m , . . . d m+n _i is weakly decreasing. Let 
d! = d m and for alH G {0, . . . , d'} define 

Xi = l + #{di,...,d m -i>d'-i}, 

V% = #{dm, ■■■ , d m+n -i > d' - i}. 

Let 

d' 

= U^- 1 - - !] x [yi-i,Vi-i + 1] U [xi - l,Xi] x [yj_i,yj]. 
i=i 

By drawing the general diagram of T(V), one immediately sees that the union above is a union of 
backwards L shapes that are piled on top of each other, and meet at the endpoints. This means 
that T(V) is a ribbon polyomino in Ribbon mjn . The bounce path is easily seen from the segments 
it is made up from, and we have 

bounce(T(:P)) = (yd'+i-i - Vd'-i, x d '+i-i - Xd'-i, • ■ • , yi - Vo, %i ~ »o)- 
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The bounce weight is 

d! 

parabounce(T(:P)) = ^ i(yd>+i-i - Vd'-i + ^d'+i-i - x d ,_i) 

i=l 

d' 

= • • • ' dm+n-i = i} + #{di, d m -i = i}) 

i=l 

= di + . . . + d m+n _i = area('P). 

So we have that T(V) G -B m ,n- Furthermore, the statistics get swapped: 

(area(T(V)), parabounce(T(7 : '))) = (parabounce(P), area(V)). 

The inverse of T is straightforward to give: Suppose the V G -B m ,n with bounce('P) = (ci, C2, 04, ... , C2 m +i) 
Then T _1 (P) is the unique polyomino "P' G ^4 m ,n with 

diaglen(7>') = (1 C2 , 2 C4 , . . . ,m C2m , (m + l) C2m+1 , m^" 1 " 1 , . . . , l Cl ). 

We omit the remainder of the details. □ 



Example 4.8. An example of the bijection from Theorem ^.7. In this case m = n = 6. The se- 
quence of diagonal lengths (di, . . . , d m+n _i) = (1, 2, 3, 3, 3, 4, 3, 3, 3, 2, 1) gi?;e (xo, . . . , 24) = (1, 1,4, 5, 6) 
and (yo, • • • , 2/4) = (0, 1, 4, 5, 6). Fill the cells according to T to get a ribbon polyomino, T(V), whose 
bounce path steps are weighted by the diagonal lengths ofV. 

T(V) = [0,0] x [0,1] U [0,l]x[0,l] U [0,3] x [1,2] U [3,4] x [1,4] 

U [3,4] x [4,5] U [4,5] x [4,5] U [4,5] x [5,6] U [5,6] x [5,6]. 
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4.3. The sandpile model on K n , special parallelogram polyominoes, and Haglund's 



bounce statistic. In this subsection we will show a connection (Theorem 4.13) between a class 
of polyominoes, the sandpile model on the complete graph K n having one sink, and Haglund's 
bistatistic (area, hagbounce) on the set of all Dyck paths. We first need to introduce some notation 
relevant to Haglund's statistics. 

A Dyck path of semi-length n is a path from (n, n) to (0, 0) that does not go above the main 
diagonal and takes steps in {s = (0, — l),w = (—1,0)}. Let Dyck n be the set of all such Dyck 
paths of semi-length n. A general D G Dyck n may be represented as a word D = d\d2 ■ ■ ■ d n where 
di G {s, w}. Given D G Dyck n , let area(D) be the number of complete unit squares contained 
between D and the diagonal line x = y. (The shaded triangular regions adjacent to the diagonal 
are not counted.) 

Haglund's bounce path for a Dyck path D is the path from (n, n) to the origin, starting with an 
initial south step, and turning from south to west, and vice versa every time the path meets the 



Dyck path D or the line x = y. If Haglund's bounce path on D is s' 



tt(i) 



w »(l) s a(2) w a(2) . 
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then we will write HagBounce(D) = (a(l), a(2), . . . , a(k)). Haglund's bounce statistic, in this setup, 
is hagbounce(D) = a(2) + 2a(3) + . . . + (k - l)a(fc). Let 

H n (q,t)= q 3re ^ D h hagbounceiD \ 

DeDyck„ 

Example 4.9. Let D = sswswswwsw = s 2 w 1 s 1 w 1 s 1 w 2 s 1 w 1 G Dyck 5 . Then area(D) = 3, HagBounce(D) 
(2, 2, 1), and hagbounce = 2 + 2(1) = 4. 




Recurrent configurations of the sandpile model on the graph K n were studied in Cori & Rossin [8] 
and are classified in terms of parking functions. (See also Knuth |19j). Let D n be the directed graph 
on the vertices V = {vq, . . . , v n -i} with single edges E = w j) : i > and < j < n — 1}. (In 
other words vertex vq is a sink.) The set of all stable configurations is 

Stable(D n ) = {(xi, . . . ,x„_i) : < x.j < n — 2 for all 1 < i < n — 1}. 

An integer sequence (ti, . . . , i n ) is a parking function if there exists a permutation 7r of {1, . . . , n} 
such that i > t^m for all 1 < i < n. Cori and Rossin [8, Prop. 2.8] proved that a stable configuration 
(xi, . . . , x n -i) is recurrent iff (n — 1 — x%, . . . , n — 1 — x n -i) is a parking function. Let Park n be the 
set of all parking functions of length n. The recurrent configurations are 

Rec(D n ) = {x G Stable(Z) n ) : (n — 1 — xx, . . . , n — 1 — x n _i) G Park n _i}. 

Let us call a configuration x G Stable(Z) n ) sorted if it is weakly decreasing, and let Rec^(D n ) = 
{x G Rec(D n) : x is sorted}. For example, 

Stable(D 3 ) = {(0,0), (0,1), (1,0), (1,1)}, 
Rec(D 3 ) = {(0,1), (1,0), (1,1)} 
Rect( J D 3 ) = {(1,0), (1,1)}. 

Using Cori & Rossin's [8J classification we have: 

Rect(D n ) = {(xi,... ) : Xj > n - 1 - i for alH G {1, . . . ,n - 1} and x% > • ■ ■ > x n _i}. 

Section 2 defined a canonical toppling process for every recurrent configuration. This canonical 
toppling process was an ordered set partition of the vertices and recorded the order in which vertices 
toppled in parallel. We now extend the same definition to recurrent configurations of K n : given 
x G Rec(if n ), let CanonTop(x) = (Qi, Q2, ■ • • , Qk) be the ordered set partition of {1, . . . , n — 1}, 
whereby vertices in Qi topples at time i. 

Example 4.10. Let x = (4,4,3,2,0) G Rec^i^)- The canonical toppling process of this configu- 
ration happens as follows: 

x + (1, 1, 1, 1, 1) = (5, 5, 4, 3, 1) -> (1, 1, 6, 5, 3) -> (3, 3, 2, 1, 5) -»■ (4, 4, 3, 2, 0) = x. 

A dot above the number denotes an unstable vertex that will topple. Thus CanonTop(x) = ({1, 2}, {3, 4}, {5}). 

Given x G Stab\e(K n ), define the diagram of x as 

n-l 

diag(x) = [n — 1 — i, 2 + X{] x [n—l — i,n — i] Cl 2 . 

i=i 
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diag(x) 




Figure 4. On the left is shown P = diag(x) where x = (5, 5, 3, 2, 2, 1) G Rec^(K 7 ). 
The Dyck path dyck(P) = s 2 w 2 sws 2 wsw 2 is shown on the right. 

Diagrams of sorted recurrent configurations in K n are precisely those diagrams which are poly- 
ominoes in Para nj „,_i. Let 

LowerPara n n _i = {diag(a;) : x G Rec T (K n )}. 

Every element P G LowerPara n)ra _i is uniquely described by a Dyck path dyck(P) of semi-length 
re — 1 which is the path of the boundary of P from (re, n — 1) to (1, 0). See Figure[4]for an example 
of this simple transformation. 

Theorem 4.11. Let u G Rec\K n ), P = diag(re) G LowerPara n) „_i and D = dyck(P) G Dyck n _ 1 . 

Then the following are equivalent: 

(i) CanonTop(u) = {Q x , Q k ), 

(ii) Bounce(P) = {\Q\\, |Qi|, . . . , \Qk\,\Qk\), and 

(iii) HagBounce(D) = (|Q X |, \Q 2 \, . . . , \Q k \). 

Furthermore: 

(iv) area(P) = (a?i + . . . + x n ) - (re - l)(n - 6)/2 = area(D) + 2(re - 1). 

Proof. The only non-trivial proof is showing the equivalence of (ii) and (iii). The Haglund bounce 
path of a Dyck path has a slightly different definition to the bounce path in a parallelogram 
polyomino. However, if one shifts the Dyck path one unit to the right, then the definitions coincide 
for this particular class thereby giving equivalence. This equivalence is illustrated in the following 
diagram: 



V 



D = dyck(P) 




The arrowed path in the diagram on the left shows that Bounce(P) = (2, 2, 2, 2, 1, 1). The arrowed 
path in the diagram on the right is Haglunds bounce path (as it is defined on Dyck paths), and so 
HagBounce(L>) = (2,2,1). □ 

Corollary 4.12. Let u G Rec^(K n ), P = diag(u) G LowerPara„, ;n _i and D = dyck(P) G Dyck n _ 1 . 
Then parabounce(P) = 2(hagbounce(P) + re — 1). 



Define the polynomial 



S n {q,t)= gareaCP^parabounceCP) 
PGLowerPara,, „_i 



Theorem 4.13. S n (q,t) = (qt) 2 ^- 1 ^ H n ^(q, t 2 ). 
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Proof. We have 



S n {q,t) = Yl g area(P) t parabounce(P)) 

PeLowerParan.n-i 

^area(D)+2(n-l)^2(hagbounce(D)+n-l) 

DeDyck,^ 

= (qt) 2 ^H n ^(q,t 2 ). 



□ 



One is led to Conjecture 4.1 by considering two configuration statistics pertaining to (sorted) 
recurrent configurations on K n . For the sandpile model on D n , the area of a Dyck path D is related 
to the number of grains in the mapped sorted recurrent configuration u since 



n-l 



Ui. 



i=0 



i=l 



Also, the Haglund bounce statistic of a Dyck path D may be interpreted as the sum over all vertices 
of the number of times a vertex is observed stable before it topples during the parallel toppling 
process CanonTop(-u). More precisely, all vertices in Qi are observed i times as stable: during the 
initial toppling of the sink and the toppling of each set Qi, . . . Q%~\. Conjecture 4.1 is the result of 
translating these configuration statistics to the sandpile model on the graph D m ,n- 



Appendix A. Appendix 

A.l. Two characterisations of parallelogram polyominoes. The bounce path may be used 
to give the following characterisation of parallelogram polyominoes. 

Characterization A.l. Let < a\ < • • ■ a m —\ < n and < b\ < ■ ■ ■ < b n < m be two weakly 
increasing sequences of integers. Let P be the intersection of the Young diagram whose corner is 
at (m, 0) having heights (n, 1 + a m _i, . . . , 1 + ai) from right to left, and the Young diagram whose 
corner is at (0, n) having widths (1 + b n , . . . , 1 + b\) from top to bottom. Then P G Para^^ iff there 
exist sequences (xo, ■ ■ ■ , Xfc+i) and (yo, . . . , yk+i) with 

(i) = xo < x\ < ■ • • < Xk = m — 1 < Xk+i = m, 

(ii) = y < yi < y 2 < • • • < yk+i = n, 

(hi) 1 + dj G {1 + Uj+i, ■ ■ ■ j Vj+2} for all 1 + Xj < i < Xj + \ and < j < k, 

(iv) 1 + hi G {1 + Xj, . . . , Xj + i} for all 1 + yj < i < yj+\ and < j < k. 

Example A. 2. The following polyomino in Paragj is formed from the sequences (a\, . . . ,as) = 
(2, 3, 4, 4, 5, 6, 6, 6) and {b\, . . . , 67) = (0, 0, 2, 4, 4, 6, 8). From this one has the sequences: (xo, . . . , £5) = 
(0, 1, 4, 5, 8, 9) and (yo, . . . , 2/5) = (0, 2, 3, 5, 6, 7) as illustrated in the diagram: 
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2/5 = 7 

2/4 = 6 > 

2/3 = 5 

2/2 = 3 > 

2/1 = 2 > 
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The next polyomino is Vi from Example 1.2, and is formed from the sequences (ax, . . . , ag) = 
(0, 0, 1, 2, 2, 2) and (pi, 62, ^3) = (3, 3, 6). The x and y sequences are (xo, x\, X2, £3) = (0, 3, 6, 7) and 
(2/0,2/1,2/2) = (0,2,3). 



2/3 

2/2 



2/0 = 2/1 = >J$- 



CO b- 



o 



1 — I (M CO 



We may give an alternative non-graphical characterisation in terms of properties of the integer 
sequence instead. 

Let m and re be fixed integers. Given two weakly decreasing sequences A = (Ai, . . . , A m ) and 
[A = . . . , fi n ) let A mi „(A, n) be the collection ofmxn cells in the plane, with Aj cells missing 
from the top of column i, and jij cells missing from the right of row j. We are interested in pairs 
of partitions A and which do not 'touch' in A mj „(A,^), as illustrated in the following diagram: 



A 



m 



Example A.3. If A = (3, 2, 2, 1, 1, 0, 0) and \i = (4, 3, 1, 0) then A 7i4 (A, fi) is V\ in Example\TJ\ If 
A = (2, 2, 2, 1, 1, 1, 0) and \i = (6, 3, 0, 0) then A^ ^X, fi) is V3 in Example \l.l\ 

Given A = (Ai, . . . , A m ) and \x = (m, . . . , /x n ) contained in the m x n rectangle as mentioned 
above. Define 

A'j = #{1 < k < n : Hk > m — 1} 
H'j = #{1 <k<m: X k >n-j}. 

for all i G {1, . . . , m} and j G {1, . . . , re}. Define A' = (A' 1; . . . , \' m ) and // = (/i^, . . . , fj! n ). Notice 
that A^ is the number of missing cells between the bottom of the mxn bounding rectangle and the 
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lowest cell in column i of A mn ((0, . . . , 0),/i). Similarly fi'- is the number of missing cells between 
the right border of the bounding rectangle and the rightmost cell in row j of A mjn (A, (0, . . . , 0)). 



Theorem A. 4. Let A = (Ai, . . . , A m ) and \i = (/ii, . . . ,/J, n ) be two weakly decreasing sequences of 
non-negative integers where Aj G {0, . . . , n} and n% G {0, . . . , m}. Then A mi „(A, fj,) G Para m)n if and 
only if 



(i) fJ-n-Xi < m — i for all 1 < i < m, 

(ii) \m~fjLi < n — i for all 1 < i < n, and 
(hi) A m = p, n = 0. 



Proof. A simple examination of adjacent columns and rows of cells tells us that A m)Tl (A, /x) G Para min 
if and only if 

(i) n > Aj + A' i+1 for all 1 < i < m, 

(ii) m > Hi + fi' i+1 for all 1 < i < n, and 
(hi) A m = /i n = 0. 

Since A^ = #{/ii, . . . , \x n > m — i} we have that condition (i) is equivalent to (i') #{1 < k < n : 
Hk < m — i} > Aj for all 1 < i < m. Similarly, condition (ii) is equivalent to (ii') #{1 < k < m : 
Afc < n — i] > ^ for all 1 < i < n. 

Conditions (i') and (ii') may be further simplified by taking advantage of the fact that the 
sequences are weakly decreasing. They are equivalent to: (i") Hn—M < m — i for all 1 < i < m, and 
(ii") Am_ w < n — i for all 1 < i < n. □ 



It follows from Theorem A. 4 by setting Aj = n — (1 + a%) for all i G {1, . . . ,m — 1} and /x. 



m — (1 + bi) for alH G {1, . . . , n} that we have an alternative characterisation: 



Characterization A. 5. Let < a% < • • -a m -\ < n and < b\ < ■ ■ ■ < b n < m be two weakly 
increasing sequences of integers. Let P be the intersection of the Young diagram whose corner is 
at (m, 0) having heights (n, 1 + a m _i, . . . , 1 + a{) from right to left, and the Young diagram whose 
corner is at (0, n) having widths (1 + b n , . . . , 1 + b±) from top to bottom. Then P G Para min iff 

(i) i < b\ +ai for all 1 < i < m — 1, 

(ii) i < ai + b t for all 1 < i < n — 1, and 
(hi) b n = m — 1. 
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A. 2. More F*^(q,t) polynomials. Some more polynomials calculated for Section 4. One can 
check by hand that they are symmetric in q and t. 



zt A q 4 (1 + qH z z + qH 6 z - qH^z 



-,3+2 



2+3, 



-,5+4^2 



H 5 Z 2 



q 7 t 7 z 3 ) 



(1 - qtz) (1 - qH 2 z) Y\ i 1 ~ i 1 ~ < f^ z ) 

(ij)e{(l,2),(l,3),(l,4)} 



1 + q 2 t 3 z + q 6 i z z + qH^z + qH^z + qH 6 z - q't^z 6 \ 5 5 
-q 9 t 7 z 3 - q 8 t 8 z 3 - q 8 t 9 z 3 - q 9 t 8 z 3 - q ll t ll z 4 J Q _ 

(1 - qtz) (1 - q 2 t 2 z) Yl (l-q i t j z)(l-q j t i z) 

(i,j)e{(l,2),(l,3),(l,4),(l,5)} 



F e ,*(q,t;z) 



( 1 + q 3 t 3 z + q 2 t A z + q 4 t 2 z + q 2 t 5 z + q 3 t 4 z + qH 3 z 
+q>ft z _ g6^6 z 2 _ „7 + 4^2 _ 9 ^10+9^3 _ ^11+9^3 



+q 13 t 12 z 4 + q^'t"z 

.18+18^6 1 ^10+15^4 



+(? 16 t 14 z 5 
+(J 12 t 13 Z 4 



t 



17+14^5 



q 7 t 4 z 2 



2q w t v z a - q'H v z 6 + 



u t n z 4 



q 10 t w z 3 



q ll t 8 z 3 + q lo t w z 



15+10^4 



>z° + q LU t LO z q 



7 17 t 19 z 6 + ^17^5 



2g 9 t 10 z 3 + 



U t U z 4 



q 8 t ll z 3 



q 9 t ll z 3 



+q 15 t W z 5 + q 16 t 15 z 5 - q 19 t 17 z 6 - q 13 t 13 z 4 - q 21 t 21 z 7 
-q 7 t 10 z 3 - qH e z 2 - qH 7 z 2 + 2 g 1 W - q 9 t 8 z 3 - q 7 t^z 3 



+q w t l2 z 4 + q 8 t 8 z 3 



q ll t 7 z 3 + 2q LZ t L1 z q 



^ 10 t 8 z 3 



12 t 12 Z A + q l2+10 z A 



+3g 

_^17 t 18 z 6 + , 
+(? 10 t 13 z 4 + , 

+q 13 t U Z 4 - , 
_^10 t 7 z 3 



g 8 t 10 z 3 + g ll t H z 4 + 



,14 t 16 z 5 



,10^4^4 
,13^4 



Q 18 t 17 Z 6 



- g 4 t V - 

g 11^3 z 4 



q 8 t 9 z 3 



3q 9 t 9 Z 3 



ji9 t 16 z* + (7 15 t 15 z 5 + q 14 t 10 z 4 - (7 6 t 5 z 2 



q i6 t i9 z 6 
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(.W fc \ (1,5), (1,6), (2, 3) / 



Here we list the polynomials F m ^ n (q,t) for small values of m and n. We use F m ^ n {q,t) 
(qt) k (m itj ) <ij< K to abbreviate 



l (q,t)= £ ™ h] q k+ H k +l. 
0<i,j<K 
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\ 1 1 / 

\l 0/ 



F 2 ,5(q,t) = (qt) 
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D 1 1110 
110 
1 0, 
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